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Abstract 

We calculate the pion decay constant near the critical temperature of the 0{N) nonlinear 
sigma model in the large N limit. Making use of the known low temperature behavior, 
we construct a Fade approximant to obtain the behavior of fn(T) at all temperatures. 
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According to folklore a nonlinear sigma model, defined by the Lagrangian 

i^ = -\{d,7.f-\{d,a)\ (1) 

together with the constraint 

7.\x)+a\x) = fl, (2) 

does not have a phase transition at any finite temperature. In other words, due to the 
constraint the magnitude of the a field is fixed at the value of and hence the symmetry 
is broken no matter what the temperature is. However, this argument is wrong. A recent 
analysis shows explicitly that a symmetry restoring phase transition occurs at 

= 12/,V(iV + 2) (3) 

to leading plus subleading order in the expansion, where — 1 is the number of pion 
field components. 

One may ask how it is possible to restore the 0{N) symmetry when the constraint 
forces the symmetry to be broken by the ground state. The answer is that a phase 
transition takes place not because the symmetry of the ground state changes, but because 
the probability of occupying the ground state decreases as the temperature increases, 
becoming infinitesimal at a transition temperature. For a familiar example, consider a 
ferro magnet. The magnitude of the spin at each site is fixed by the constraint = 
5'(5'+l). The ground state of this system is of course one where all the spins line up 
in a fixed direction, breaking the rotational symmetry of space. However, this system 
certainly goes through a symmetry restoring phase transition as the temperature is raised. 
At high temperatures the amount of energy put into the system is simply too great to be 
accommodated by varying the phase of the ground state alone. 

Given that there is a phase transition in the nonlinear 0{N) sigma model, it is im- 
portant to ask how the pion decay constant changes as the temperature increases. One 
must carefully define what is meant by the pion decay constant at finite temperature. 
A reasonable definition is that it measures the strength of the coupling of the pion, the 
Goldstone boson, to the longitudinal part of the correlator of two axial vector currents. 
Following we take as the definition 

/,^(T) = 21ini r ^PL{qo,q=0). (4) 

Jo 

Here Pl(q') is the longitudinal part of the spectral density for the retarded response func- 
tion of two axial vector currents, where 

= Crdf^TTa - -Kad^a (5) 

in the sigma model. To give a more physical meaning, consider first what pions are 
in a finite temperature system. Pions are Goldstone bosons associated with the broken 



2 



chiral symmetry, and in a finite temperature system represent the energy stored in a 
varying phase structure of the ground state. Consequently, the pion contribution to the 
correlation function must depend on the probability of the system to be in the ground 
state. That probability is relatively large at low temperatures and so the pion decay 
constant is not much different from the zero temperature value. As the temperature goes 
up, the probability of being in the ground state diminishes, and so does the pion decay 
constant. At the phase transition temperature no pions are present in the system since 
no energy is stored in the varying phase of the ground state, and the pion decay constant 
as well as the scalar condensate becomes zero. This all assumes the transition is second 
order, of course. 

The constraint is local and difficult to deal with directly. It is advantageous to repre- 
sent it with a Lagrange multiplier and an auxiliary field X{x) as is done in P, 0]. In 
Euclidean space, used throughout this paper. 
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(d„TTy + (d^ay + m'(Tr' + a') + a(-K' + a') . (6) 



The auxiliary field A over which we must integrate to get the partition function must lack 
a zero frequency and momentum component 

j£xdTX{x) = 0, (7) 

otherwise TT'^{x)+a^{x) will be identically zero. The Lagrange multiplier together with 
A enforces the averaged constraint 

(tt^) + {a') = fl . (8) 

The next step is to obtain an effective action involving only the auxiliary field A by 
integrating out the tt and a degrees of freedom. At a first glance, this seems like a simple 
task since the above Lagrangian (^) is bilinear in tt and a. One is tempted just to carry 
out the necessary Gaussian integrals in the partition function 

Z = j [dTT][da][d\] exp (^j d^xdrC^ . (9) 

However, doing so gives a wrong result. A Gaussian integral 

J d^x exp [-XT Ax] = 7r"/2 Det-^/^(A) , (10) 

is valid only if the determinant of A is positive definite. When = 0, the action is 

not positive definite since J d^xdr C vanishes when both tt and a are constant. Hence, 

zero modes exist when = 0, and if the partition function is to be well-defined the 
integration must not be done over constant field configurations. 
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To carry out the integration over only the non-zero modes it is convenient to separate 
the constant zero modes exphcitly from the non-zero ones. Due to the 0{N) symmetry 
one can always set (tt) = and (cr) = v. Shifting a hj v then yields 
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{d^nf + {df,af + m^{TT^ + (r^ + 2va) + ^(Tr^ + + 2w) (11) 



where J d^xdr X = is used and a constant term -m^vV2 is discarded. The condensate 
V is determined by the condition 

(a) = 0. (12) 

To satisfy the condition (0), the action must not contain terms odd in a. Since A has 
no k = component, the only such term in the Lagrangian (|TT]) is —m^va. The condition 
then implies that and v cannot both be nonzero. As explained in [jl|, there are two 
possibilities.^ 

1. = and v ^ 0. This implies that the symmetry is broken and there are massless 
Goldstone modes. Since the symmetry is broken at zero temperature, this must 
correspond to low temperatures. 

2. 7^ and v = 0. This implies that the symmetry is restored, and hence there are 
no massless Goldstone modes. This must correspond to high temperatures. 

The constraint (P) and the condition (p!2D must be imposed after the partition function is 
calculated with arbitrary and v. 

Carrying out the Gaussian integrations over tt and a requires completing the square. 
In order to do so, a is shifted once more by a solution of the equation of motion 



ax{x) = 2iv J dz Gx{x, z) X{z) , (13) 

where 

Gx = l/i-d^ + -2iX). (14) 
Replacing a with a+ax gives the Lagrangian 

= -2^2 A GxX + ^TT- {d^ -m^ + 2iX) tt + {d^ - + 2iX) a. (15) 

Performing the Gaussian integration over the non-zero modes of tt and a expresses the 
partition function in terms of only the auxiliary field A. 

Z = exp [-{N/2)Tt' In (-9^ + m^)} J [dX] exp (Sx) , (16) 



^There is, of course, the possibility that both and v are zero. However, if this is true, we wouldn't 
be able to satisfy both of the conditions (ph and (|3). 
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where 

Sx = -2v^ J d'xdr A Ga A - (iV/2)Tr In (^1 + , (17) 

is the effective action for the A field. The prime on the trace indicates that the zero 
eigenvalue component is missing from the operation. Note, however, that trace-log of 
[1 + 2iX/{d'^ — m^)] is well-defined even when = since the A field has no k = 
component. Consequently the trace operation in the effective action Sx can include the 
k = component without any error. 

At finite temperature, as argued in the residue of the pion pole in the retarded 
response function for the axial vector current is proportional to the square of the pion 
decay constant. Hence, calculating the spectral density (see, for example, 0) 

in the zero momentum and small frequency limit yields f^{T). We now calculate the 
spectral density by taking the discontinuity of the Euclidean Green function. Then, from 
the coefficient of the massless pion pole, f^(T) is obtained. 
As before, one first shifts the sigma field twice to yield 



= [d;rr,{x)[v + a{x) + ax{x)]-n^{x)[d;a{x) + d;crx{x)] 
X {dyMy)[v + cr{y) + ax{y)] - Tr.iyWMv) + d'^^y)]] ■ (19) 



This expression is to be averaged with the effective Lagrangian given in eq. (|15|). Since 
the latter is quadratic in tt and a, Wick's theorem trivially works with the propagator 

{TTaix) TTbiy)) = {crix) (j{y))5ab = {Gxix, y))5ab ■ (20) 

Performing the tt and cr integrations is equivalent to replacing each bilinear in eq. ([T9| ) 



with Gx and yields 

{^".(x) Aliy)) = cf^dl (Ga(x, y))b,b + ( V(x, y) + B^,{x, y) + C,,{x, y))5ab , (21) 
where 

A, ,{x,y) = 2[Gx{x,y)&;,dlGx{x,y)-ff;,Gx{x,y)dlGx{x,y)], (22) 

B, Ax. y) = vld^^d^Gxix, y) ax{x) + &;^dlGx{x, y) ax{y) 

- dlGx{x,y)dlax{y)-dlGx{x,y)dlax{x)], (23) 

C^^(x, y) = dldlGx{x, y) ax{x) ax{y) + Gx{x, y) d^axix) dlax{y) 

- dlGx{x, y) d;ax{x) ax{y) - d;Gx{x, y) ax{x) dlax{y) . (24) 



Here the facts (tt) = (a) = {a\) = are used to simplify the result. The expression ( PT| ) 
can be now integrated over the remaining A field with the effective action (|1^) to yield 
the axial vector current correlation function. 

First, consider the unbroken phase where v = 0. In this phase the condensate is gone 
and no pion is present in the system. Consequently the pion decay constant must be zero. 

The expression (PT|) is exact. To use it to calculate /^(T) in the broken phase where 



f 7^ and = 0, one must now distinguish temperatures much lower than the transition 
temperature and temperatures close to Tc. 

First, consider temperatures near Tc so that f^(T) <^ NT'^. At these temperatures 
the effective action ( |l7|) is of 0{N). Hence, the propagator for the A field is of 0{1/N). 
Ignoring A field contribution all together in eq. (^) yields for the (00) component of the 
Euclidean Green function 

r d^l 

Gf^{k) = kokov'G{k) + W^^oi^lo - h) G{1) G{k - I) S^b , (25) 

where G{k) = The first term in eq. (^) clearly contains a pion pole. The second 

term remains finite as A; — after a suitable zero temperature renormalization. Conse- 
quently only the first term contributes to f^(T). 

To obtain the pion decay constant at finite temperature one first takes the disconti- 
nuity across the real ko axis and then divides the result by fcg. This results in 

PAT) = v\T) = p^-^T'- T^T^. (26) 

Here the value of f^(T) near Tc is extracted from There are 0{v'^/N) corrections to 
this result from further expansions of v'^d^d^Gx, B^^ and G^u terms in eq. (|2T|) . There is 
also an 0{T'^/N) correction from further expansion of A^iy. 

The behavior of the pion decay constant at low temperature was first obtained by 
Gasser and Leutwyler for QCD |Q and subsequently verified using a different method by 
Eletsky and Kogan 0. For the 0{N) sigma model the result isQ 

TV — 2 

f^{T) = fl - -^T' ; T « T, . (27) 

We have also been able to reproduce this result using the methods developed in this paper. 

To compute the behavior of (T) at arbitrary temperature is not much more involved 
than the calculations performed in this paper. One just needs to use the full propagator 
(pop for the A field instead of the limiting one. Such a numerical evaluation of /^(T) 
as a function of T would be in itself an interesting project. However, it is also clear 
that no surprising feature will emerge from this calculation since the A propagator is 

^The group SU{Nf) x SU{Nf) is isomorphic to 0{Nj) only for Nj ~ 2 and this hmits the quantitative 
comparison to two quark flavors. 
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Figure 1: Plot of the Pade approximant eq. ( p^ with = 4. The broken hnes are hmiting behaviors 
given by eqs. (p6|) and (p7|). 



always small compared to the temperature. Hence, for most applications that are not 
terribly dependent upon fine structure, the following Pade approximation between the 
low temperature result (|27|) and the near result (|26|) should suffice. 



1 - [4/(iV + 2)] (TVT2) (1 - TVT2 



(28) 



The above Pade approximation is plotted in Figure 1 for the case = 4. 

A natural extension of this work is to include the vector and axial vector mesons, 
p,u,ai, and recompute the temperature dependence of the pion decay constant. This 
would allow for a more direct and meaningful comparison with lattice gauge computations. 
However, lattice gauge calculations of the pion decay constant at finite temperature so far 
have extracted it using the Gell-Mann Oakes Renner relation 0. It is not clear that this 
is a good definition at finite temperature, nor whether it is the same as the definition in 
eq. (4). In addition, the lattice calculations have used the quenched approximation which, 
while of pioneering form, may be missing some essential physics of the chiral transition. 
Further studies in both theoretical and lattice calculation are clearly called for. 
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